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Abstract 

The quark recombination mechanism is proposed to explain the asymmetry be- 
tween production rates of and of D~ in their inclusive productions, and also 
asymmetries for other charmed hadrons. These asymmetries are observed in ex- 
periment and are called as leading particle effects. In this work we give a general 
analysis for contributions of quark recombination to these asymmetries. The con- 
tributions consist of a perturbative- and nonperturbative part. We perform two 
types of factorization by considering the produced hadron with large transverse 
momentum and by taking charm quark as a heavy quark, respectively. In the 
case of large transverse momentum the effect of quark recombination is the stan- 
dard twist-4 effect. We find that the contributions are parameterized with four 
nonperturbative functions, defined with four quark operators at twist-4, for initial 
hadrons without polarization. By taking charm quark as a heavy quark the fac- 
torization can be performed with the heavy quark effective theory (HQET). The 
effect of quark recombination is in general parameterized by eight nonperturba- 
tive parameters which are defined as integrals with matrix elements of four quark 
operators defined in HQET. For unpolarized hadrons in the initial state, the pa- 
rameters can effectively be reduced to four. Perturbative parts in the two types of 
factorization are calculated at tree-level. 

PACS numbers: 12.38.Bx, 12.39.Hg, 13.85.Ni 
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1. Introduction 



Because nonperturbative physics of QCD is still unclear, a prediction of inclusive hadropro- 
duction of a hadron can not be made in general. However, if the produced hadron has a large 
transverse momentum kt, one can predict the differential cross section. According to the factor- 
ization theorem in QCD (see e.g., [1]), the prediction takes a form as a convolution with parton 
distributions for partons in the initial hadrons, fragmentation function of a produced parton decay- 
ing into the produced hadron, and perturbative functions for partonic processes. The prediction is 
based on the leading terms in the so called twist expansion. In the leading terms twist-2 operators 
are used to define nonperturbative parts, i.e., parton distributions and fragmentation functions. 
It can be shown that the next-to-leading orders are suppressed by certain power of kf relatively 
to the leading order. At the next-to-leading order, i.e., at twist-4 level, various partonic processes 
can happen. The nonperturbative effects are parameterized by matrix elements of twist-4 opera- 
tors. At this order, a produced meson can be formed by recombination of two partons in a initial 
hadron. This is the so called parton recombination model [2], which is proposed to predict the 
inclusive production of a meson with low kf. The recombination can also happen with one of the 
two parton involved in a partonic process. This has been studied in [3] . 

If we take the produced hadron containing one heavy quark in a hadron collision, for 
example, D"*", and neglect possible c-quark content in the initial hadrons, then at leading order 
of perturbative QCD the production can be understood as that a cc pair is produced via 
partonic processes like gg cc or qq — > cc and then the c-quark fragments into the meson. 
It is clear that the production rate for and for D~ is the same at this order. Experimentally 
significant asymmetries between these production rates for a class of c-flavored hadrons have been 
observed and have been called as leading particle effects. These effects are observed in fixed-target 
hadroproduction experiments for the production of c-fiavorcd hadrons [4,?,?,?,?] and also in pho- 
toproduction experiments [9-11]. Significant asymmetries are observed for — , D'^ — D^^ and 
— while for Df and A+ the asymmetries with large errors are consistent with zero. These 
asymmetries are large when kt is small, for and D~ it can be ~ 0.7 in the forward direction 
in hadroproduction. Theoretically, various explanations exist [12,13], e.g., these asymmetries can 
be generated by intrinsic charm in the initial hadrons, but the predicted asymmetries are smaller 
than those in experiment. It should be noted that these asymmetries can also be generated at 
next-to-leading order of perturbative QCD [14-17], the predicted asymmetries are too small to 
explain the observed. 

Recently, it has been suggested that these leading particle effects can be explained with quark 
recombination [18-20]. In E791 experiment [4] a vr" beam is scattered on a nuclear target, a 
D~ meson can be produced through quark recombination of a produced c-quark with the valence 
quark d in tt" after hard parton scattering, while a D'^ meson can be produced through quark 
recombination of a produced c-quark with the sea quark d in tt". Because the distribution of d- 
and d quark in 7r~ are significantly different, it results in that more D~ are produced than 
in the forward direction. Motivated by NRQCD factorization [21], the authors of [18-20] take a 
factorization approach for quark recombination, in which a c-quark combines with a d quark in 
the case for meson and the combination cd is in ^Sq or ^Si state, i.e., the pair cd is formed 
as a pseudoscalar state or vector state, the cd is then inclusively transformed into the D'^. If 
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the d quark was a heavy quark, the above states would be dominant in the quark recombination. 
However, the d quark is a hght quark, for it a NRQCD description is not correct. It is possible 
that the cd is produced in a scalar, pscudovector and tensor state and then transformed into D^. 

If one takes c-quark as a heavy quark, a factorization different than that mentioned at the 
beginning can be performed with the Heavy Quark Effective Theory (HQET) [22]. In the heavy 
quark limit mc — > oo, the c-flavored hadron produced in a collision will carry the most momentum 
of the produced c-quark which is transformed into the hadron after its production. This suggests a 
factorization for the production of the c-fiavored hadron, and in this factorization the production 
rate can be written as a product of the production rate of c-quark and a matrix element defined 
in HQET, corrections to this are suppressed by certain power of m~^. This factorization was 
used to predict heavy quark fragmentation functions [23] and spin alignment of heavy meson in 
its inclusive production [24]. With this factorization predictions for these two cases are in good 
agreement with experiment at e+e~ colliders [23,24]. 

In this work we will take this factorization to analyze the possible source for generating leading 
particle effects. We find that in general the cd pair can be first formed into the states of scalar, 
pseudoscalar, vector and pseudovector, from these states a D'^ can be generated by emitting 
unobserved states. Our final results consist of 8 nonperturbative parameters, which are defined 
with matrix elements of HQET. Corrections to our results are suppressed by the inverse of heavy 
quark mass. These parameters can be calculated with nonperturbative methods. For initial 
hadrons without polarizations, these 8 parameters can be effectively reduced to 4 because some 
perturbative coefficients of the parameters in their contributions to the production rate are same. 
If the initial hadrons are polarized, all perturbative coefficients can be different, the number 
of parameters can not be reduced. In the approach of [18-20] the contributions from quark 
recombination for unpolarized initial state also contain 4 nonperturbative parameters, but their 
interpretation is different than our 4 parameters obtained from the effective reduction of the 8 
parameters. Our results hold not only for low kt, but also for high kt- If k^ is enough large so that 
all quark masses can be neglected, the standard factorization mentioned at the beginning can be 
performed. In this work we also perform an analysis for this case. This is similar to the analysis 
of twist-4 effects in deeply inelastic scattering [25], however, our task here is simple, because 
we only need to analyze the contributions from twist-4 operators, which are defined with 4-quark 
operators. These twist-4 operators are corresponding to those 4-quark operators in deeply inelastic 
scattering. In our final results for large kt, there are four nonperturbative functions defined with 
4-quark operators separating along the light-cone. It should be noted that it is rather difficult to 
perform a factorization for the relevant processes, in which hadrons are in the initial- and final 
state. We will assume that the factorization holds and especially that the parton model can be 
used for initial hadrons. With this assumption we need only to consider inclusive production of a 
hadron in a parton collision and focus on the final hadron. 

Our work is organized as the following: In Sect. 2. we give our notations and analyze the 
quark recombination in the case with large kt where the kt is so large that all quark masses can 
be neglected. The contribution to the differential cross-section is a convolution of perturbative 
functions for partonic process with nonperturbative functions defined with four quark operators. In 
Sect. 3. we analyze the quark recombination with the mentioned factorization with HQET, in which 
nonperturbative parts are defined as matrix elements of HQET. There are eight nonperturbative 
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parameters characterizing different states of a quark pair which is transformed into the produced 
hadron. Wc also show in detail that these 8 parameters in the case with unpolarizcd beams can 
be effectively reduced to 4. Discussions of our results are given. Sect. 4 is our summary. 

2. Quark Recombination at Large kt and Twist-4 Effect 

We will denote a heavy quark as Q and a heavy meson containing one heavy quark Q as 
Hq. The heavy meson Hq has a valence quark q. As discussed in the introduction we will 
consider the production of Hq in a parton collision and a initial parton is a valence parton of Hq. 
Neglecting possible heavy quark content in initial hadrons, the relevant parton process for quark 
recombination is: 

9{pi) + q{p2)^HQ{k)+X, (1) 

where momenta are given in the brackets. We can always divide the unobserved state into a 
nonperturbatively produced part and a perturbatively produced part Xp, i.e., X = X^ + 
Xp. At leading order of ag, Xp is just a antiquark Q. The scattering amplitude for the quark 
recombination can be written: 

T = l ^^Aij{p^,p,,k,,p,) I d^Xie-*-i(i/Q + Xjv|g,(xi)g,(0)|0), (2) 

where the indices i,j stand for color- and spin indices. Q{xi) and q{0) arc the fields of Q and 
q respectively. Aij is the amplitude for g{pi) + q{p2) — ^ Q*{ki) + q*{K) + QiPs): * means that 
the states are not on-shell. If the states are on-shell, then the amplitude is u{ki)iVj{kq)Aij with 
on-shell conditions for corresponding partons. The Feynman diagrams for the amplitude A^j are 
given in Fig.l. For the process we define the variables: 

s = (pi+P2)', i^{p2-kf. (3) 




(d) (e) 



Fig.l. Feynman diagrams for the amplitude A. 



The contribution from the process in Eq.(l) to the differential cross section is: 
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d^k d^k 



where we use the subscribe R to denote the contribution from quark recombination, the average 
of spin and color of initial partons is implied. If Hq has nonzero spin, the summation over the 
spin is understood. Using translational covariancc one can eliminate the sum over Xj^. We define 
a'ff^ as the creation operator for Hq and we obtain: 

1 d^k f d^ps f d^ki d^k2 ^ o.t^ u \ o\ 

d''x,d^X2d''xze-'^'--''^'^^--^-'^'-^'mk{m^ (5) 
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with kz — pi + p2 — kl — pz s& the momentum of q after the hard scattering 
in the amplitude Aij{pi,p2iki,pz}- This contribution can be conveniently represented 
by the diagram in Fig.2., where the lower part represents the perturbative part, i.e., 
^ii(PiiP2i fciiP3)(7°^^(PiiP2! ^2iP3)7°)fc«) the upper part with the black box represents the non- 
perturbative part, i.e., the Fourier transformed matrix element in the second line of the above 
equation. 




Fig.2 

Fig.2. Diagram representation for Eq.(5). In the figure k^ — ki — k2-\- kz, the dash line is the cut, 
the thick lines below the black box are for heavy quarks. 



The factorization can be performed by decomposing the color- and Dirac indices of the matrix 
element. The decomposition of color indices of the matrix element can be performed in a straight- 
forward way and one can identify that the Qq pair is in a color-singlet or a color-octet. We denote 
the color-singlet- and color-octet contribution as da^^^ and da^^^ respectively and 
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ddR^da^^^ ^daf. (6) 

The contributions read: 

(1) 1 d^k r d^ps r d^ki d'^k2 
" 2i (27r)32P J (27r)32p0 ' J (2^ (2^ 

1 /■ 

•gTrAj(Pi,P2, A;i,P3)Tr(7°^^(Pi,P2, A;2,P3)7°)w j d'^Xid'^X2d'^x^ 

.g-ifei..,+ifc....-i.3-*:3^0|g,(0)gK^2)aLjA;)aHjfc)gi(xi)(Z,(a;3)|0), (7) 



(8) 



d^k f d^ps f d^ki d'^k 



^ ^ 2s (27r)32A;0 J (27r)32p0 J (27r)4 (27r) 

~IVr%(pi,p2,A;i,p3)Tr(rV^^(pi,P2,A;2,P3)7°)M / d^xid^X2d^xs 
•e-''=i-"i+''=^-"^-'"^-'=^(0|gfc(0)r"QK^2)aL^(A;)aj^^(A;)gi(xi)r%(x3)|0). (8) 

In the above equations the traces are taken over color indices and the indices i,j, k, I stand only 
for Dirac indices. 

In this section we will consider the case where the transverse momentum kf of the produced 
hadron is very large, and all quark masses can be neglected. To evaluate the contributions we take 
a frame in which the ^-direction is the moving direction of Hq and define a light-cone coordinate 
system in which a vector V is expressed with components V^{ii = 0, 1,2,3) in usual coordinate 
system as: 

= {v\ V-, v\ V-) = [^-^, v\ j . (9) 

We also introduce two light-cone vectors: 

= (0,1,0,0), r = (1,0,0,0). (10) 

The momentum k reads: 

Ml 

k^ - (^+, k-, 0, 0) - (^+, 0, 0) « (^+, 0, 0, 0) = k+F. (11) 

The x-dependence of the matrix elements in Eq.(7) and Eq.(8) is controlled by different scales. 
The x~ dependence is controlled by the large scale k'^, while the x+-dependence is controlled by 
the small scale k~ or Aqqd, the x^— and x^-dependence is controlled by Aqcd- Therefore, as 
an approximation one can neglect the dependence of x~ , x^ and x^. Corrections to this can be 
systematically added. This results in the so called coUinear expansion. The leading order of the 
expansion is to set 

kf ^ {zik-^, 0,0,0) ^ Zik+F, fori = 1,2, 3, (12) 
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in the amplitude A and A^. Taking the color-singlet in Eq.(7) as an example, we obtain: 

(1) 1 d^k r d'^ps r dzi dz2 ,,+^2(r, \3r3 



{^\qk{^)Qi{x2n)a^HQ{k)<^HQ{k)Qi{x^n)qj{x^n)\Qi), (13) 

where = 5{q^)5{q^)5{q^). With Eq.(12) the variable z~^{i = 1,2,3) is just the momentum 

fraction of partons carried by Hq. Because the momentum conservation they are bounded, for 
example, > 1 for i = 1, 2, 3. Now the decomposition of the Dirac indices can be performed. We 
will only keep the leading terms, i.e., the twist-4 terms. The decomposition read: 

= {l~)ji{l~)ikTi^\zi,Z2,z^) + (7"75)j-i(7-75)ifeT2^^^(zi,2;2,2:3) 

+ {ta-^),,{ia-niai'\zu Z2, ^3) + ■ ■ ■ , (14) 

where ■ ■ ■ denotes the terms with twist higher than 4. The sum over the index m runs from 1 to 
2. The functions %{% = 1, 2, 3) are defined as: 

tI'\z,,Z2,z,) = — [ ^^^e-'=^(^^-^-^^-^"+^3.3-) 

1 y 1, 2, 6) J 2^ 2^ 27r 

• (0|g(0)7+(5(x^n)a]:^^a//Qg(a;^n)7+g(a;3"n)|0), 

2 V i> 2, 6) J 27r 27r 27r 

• (0|g(0)7+75Q(a;^n)ai^^a/^QQ(a;^n)7+75g(a;3n)|0), 

r3('^(^i,^2,^3) = — / ^^^e-'=^(^^-r-^^-^+^3.3-) 

3 V i> 2, 3; J 27r 2tt 2tt 

■ (0|g(0)icT+*Q(a;2n)al a// O(a;]"n)i(T+*g(a;Jn)|0), (15) 



these three functions represent the nonperturbative efi^ects related to the produced hadron. They 
are invariant under a Lorentz boost along the moving direction of Hq. The dimension of these 
functions are 2 in mass, hence they are proportional to ^cd- With these results we obtain 



, (1) i «fc f d P3 fdzidZ2, ,3 3, , 
^ 2i (2^ i (2^ ■ J — — ^'"^ '^-^^^ +P2 - ^1 - P3) 

;{Tr(7 • kiA{pi,p2,ki,p3))TY{'y ■ k2'y°A\pi,p2,k2,P3h^) ■ 1i^\zi, Z2, Z3) 



9k+ 

+IV(7 • A;i75^(pi,P2, A:i,P3))TV(7 • A;2757°^^(Pi,P2, A;2,P3)7°) • ^2 ^1, ^2, ^sj 
+Tr(ia^'%^A(p,,p2, ki,p3))Tr{ia^%^^''A''{p,,p2, A;2,P3)7°) ' '^a^'H^i, ^2, ^3)}, (16) 
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for the color-octet contributions we obtain: 
1 d^k r d^P: 



, (8) -L "-/i- r a-^p^ fdzidZ2, .3 3. . 



2A;+ 



{Tr(7 ■ A;iTM(pi,p2, A;i,P3))Tr(7 ■ k^^^T^ A^p^^p^, A;2,P3)7°) ■ ^/'^(^i, ^2, 



^3 



+Tr(7 • k^-f,T^A{p^,p2, k^,p^))^T{^ ■ k^i^i'T'^A^p^.p^, k^.p.h') ■ tI'\z,, 



+Tr(ia^%^T-A(p,,p-,, A;i,p3))IY(ia''^A;2^7°rMt(pi,p2, A;2,P3)7°) ' ^2, ^3)}, (17) 

with ki = Zik~^l. The definitions for the functions T^^^\zi, Z2, Z3){i = 1,2,3) can be obtained by 
replacing the color matrix 1 1 with T" T" in the definitions for T^^\zi, Z2, Z3){i = 1,2,3), 
respectively. 

The integration over phase space can be easily performed in the center-of-mass frame. Consider 
the integral 



{2n)^2k^J (27r)32p§ 



(18) 



because the only nonzero component of ki is kf, ki can be set zero in the 5- function. The 
integration over p^ can be done with the 5- function. We insert in the above integral the identity: 



j dz3k+d{pt +P+- zik+ - Z3k+) = 1 = / dz3V2\k\S{^ - V2{zi + Z3 



(19) 



and obtain the phase space integral as 
d^k r d^ps 



r ark. r 
J (27r)32A;0 J 



(27r)32A;0y (27r)32p0 



Q-{2ny5t{pi+p2-ki-p3) 



k+ 4(27r)2s 



dt / dzs 



{zi + ZsY 



(20) 



With this we can express our results for da/di. The remainder of the calculations are straight- 
forward. Without detailed calculation one can find that the coefficient of T^^^ is zero, the reason 
is that we consider unpolarized hadron beams here and the polarizations of initial partons are 
averaged. Due to helicity conservation, the coefficient must be zero. If we consider polarized 
hadron beams, T^^^ can lead to a nonzero contribution. With the structure of T^^^ and T^*-*^ one 
can also find that the coefficients of 7^*^*^ and T^^^ are same. Finally we have: 



dcR\^ '\ 27ra'^l f dzi dz2 dz-^ 
(s, t) 



dt 243s2 

S2 



1 



Zl Z2 Z-i [Zi + 2:3)2 (zi - Z2 + Z3){S + t) 



1^(9^1 + Z3){zi + 8z2 + ^3) + ^ [27^1^ + 297^1^2 + 29^1^3 + 25^2-^3 + 2zl 



S2 



n 4t r 

81^3(^3 - Z2) + Zi{25z2 + 81;23)J - -[8lzl + 9z3{8z2 + Z3) + Zi{90z3 - 86^2 



-144z^ + 469^1^2 - 179^1^3 - 109^2^3 - 35z^}, 



(21) 
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da 



(8) 
R 



dt 



{s,i) 



iral f dzi dz2 dz^ 



I 



(8)1 



Zi Z2 Z3 {Zi + ZsY [zi -Z2 + Z3) 



S + t 



{^[GSzI + Z3{Z2 + 62^3) + ^1(9^2 + 125^3) 



■s 

? 

s 
i 

7? 



+^ [207zl + 99ziZ2 + 331^1^3 + 83^2^3 + 124z 

At 



+— [45^3(^3 - Z2) + Zi{Alz2 + 4:5zi 
-18zl + 5^1(79^2 - 47^3) + -23(199^2 - 217^3)}. 
Similarly for the process 

^{pi) + q{p2)^HQ{k)+X, 

we have 



81^^ + 9^3(^2 + 8Z3) + zi{153z3 - U8Z2) 



(22) 



da 



R 



dt 



{s,t) 



3P 



dzi dz2 dzs 



J c(-)e|{r,^'^ + r2^n 

Zi Z2 Zs {Zi + 2:3)2 {Zi - Z2 + Z3) 



(23) 



(24) 



{ 



[i^ZiZ2 + s'^Z3{zi - Z2 + Z3) 



p (s + i) 

where eq^Cq) is the charge of the heavy(light) quark and 



_ 256 

eq 27 



^ ~9- 



(25) 



With these results one can predict the leading particle effect in hadron collision and pho- 
toproduction. We denote the anti-particle of Hq as Hq. Because of the symmetry of charge 
conjugation, the differential cross section for g{pi) + ^'(^2) Hq{k) + X is the same as given 
above. The leading particle effect is generated by the asymmetry between the distributions of q 
and q. We consider the production in the hadron collision: 



A{Pi)+B{P2) 

A(Pi) + B{P2) 



Hq{k)+X, 
Hq{k)+X, 



(26) 



where A and B are the initial hadrons, whose spin is not observed. We define the variables: 

s^(P, + P2f, t^{P2-k)\ u^{P^-kf. (27) 
The leading particle effect can be predicted as: 

}) = j dXidX2fg/A{xi) [fq/B{x2) - fq/B{x2) 
+ j dXidX2fg/B{x2) [fq/A{Xl) - fq/A{Xl) 
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^(Hq)-^{Hq) 



da 



R 



dt 
dan 

di 



S — X1X2S, t — X2t) 

s — X1X2S, t — xiu), (28) 



where dan is the sum of contributions given in Eq.(21) and Eq.(22), fa/A{x) is the distribution for 
parton a in hadron A. Similarly, one can also obtain the leading particle effect in photoproduction, 
replacing A with a photon in Eq.(26) we obtain: 



-^W - -{Hq 



J dXi [fq/sixi] - fq/sixi)] " "^(^ = ^ = ^1^' (29) 

where dan = da^j^^ + da^^\ If A is the anti-particle of B, then the leading particle effect vanishes 
because of the charge conjugation symmetry. Since the nonperturbative functions 7]^'^^ and T^^^ 
for i = 1,2 are unknown, numerical predictions can not be made. These functions may be studied 
with nonperturbative methods or extracted from experimental results, they are universal like 
parton distributions. Once they are known, one can make numerical predictions. Our results here 
hold for the case with large transverse momentum of the produced hadron, i.e. For s — > oo and 
i — > oo, but with the ratio t/s being fixed, one obtains the power behavior 

(30) 



-(Hq) ~ -(Hq) ~ (31) 



in comparison with the production at twist-2 level 

-(Hn) ~ -(Hq) ~ - 

dV ^' dV ^' s2 

hence, the leading particle effect through quark recombination is suppressed in comparison of the 
production rate by a factor s~^, which is typical for twist-4 effects. 

In [3] contributions from quark recombination to inclusive tt production is studied, partonic 
processes like gq Tig and gq — > qn are calculated by taking a wave function at leading twist for 
TT. This is different than our approach. If one only keeps the contribution from T^''^-' in Eq.(21) 
and neglects all other contributions, and then for T^*-^-* one uses the approximation of vacuum 
saturation to write it in term of the n wave function at leading twist, our approach is equivalent 
to that in [3]. It is clear that there are no reasons to neglect other contributions and to use 
the approximation of vacuum saturation. In [3] the results are given by taking asymptotic form 
of the wave function, where the integration over the momentum fraction as the variable of the 
wave function is performed analytically. This prevents us from a direct comparison with our 
results. From our results it is clear that the effect of quark recombination is characterized by 2 
nonperturbative functions, effectively for unpolarized beams, while for polarized beams there are 
in general 6 nonperturbative functions. 

3. HQET Factorization for Qucirk Recombination 

In this section we will use nonrelativistic normalization for Hq and heavy quarks Q. If the 
pole mass mg of the heavy quark Q is heavy enough, one can expand the heavy quark field Q{x) 
with the field in HQET [22]: 

Q{x)^e-^'^'^^-^-{h{x)+0{mQ')} + ---, 

Q{x) = e+'"*'3'"^ • {h{x) + C»(m^^) } + •••, (32) 
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where v is the four velocity of Hq, i.e., v — k/Mff^. h{x) is the field of HQET and depends on v 
implicitly, the • • • denotes the part of anti-quark. The fields have the following property: 



V ■ jh{x) = h{x), h{x)v ■ 7 = h{x). (33) 

For ttlq — > oo the most momentum of a produced Q is carried by Hq, where the produced heavy 
quark Q combines with other light quarks and gluons to form the hadron Hq. Using this fact 
and HQET one can expand the cross-section for production of Hq in I^qcd/^q- Following the 
proposed factorization in [23,24], the production rate at leading order of ^qcdI'^iq is a product 
of the production rate of Q with a matrix element defined in HQET, the momentum of Hq is 
approximated by the momentum of Q and Muq ~ vtiq. For example, the production of Hq via 
gluon fusion, i.e., g + g ^ Hq + X, the differential cross-section can be written as: 

da{g + g^HQ + X)= da{g + g ^ Q + X) ■ ^(0|Tr/ial, a^^/i|0). • (l + 0{^^)\ , (34) 

D TUq J 

where matrix element is defined in the rest frame of Hq. The above equation has a close corre- 
spondence to inclusive decay of Hq with the factorization, the corresponding matrix element for 
inclusive decay of Hq equals 1 because fiavor conservation of HQET, it leads to the well known 
result [22]: 

T{Hq^X) = T{Q^X)-[i + 0{^^^^)\. (35) 

At the next-to-leading order for the decay width the nonperturbative effect is represented by 
dimension 5 operators, corresponding operators for production can also be found. These oper- 
ators are bilinear in quark fields. Going beyond this order, i.e., at the order of rri^ , one will 
encounter four quark operators, these 4-quark operators represent the nonperturbative effect for 
quark recombination. 

For the contributions from quark recombination, we neglect higher orders in Eq.(32). The 
color-singlet contribution reads: 

• j (l'^xid'^X2d^x^e~^'^^'^^~^^^^'^^~^^'^'^'^~^"'''^^'^^~^^"^'^^'^^ 

•(0|^fe(0)/iK^2)alj^(A;)a^^Q(A;)^i(xi)g,(x3)|0). (36) 

Now the space-time dependence of the matrix element is controlled by the small scale Aqcd, 
refiecting the fact that the most momentum of Q is carried by the hadron Hq, the difference 
between the momentum of Q and that of Hq is order of Aqcd- At first look, one may neglect the 
space-time dependence as an approximation, used in deriving Eq.(34) or Eq.(35). This implies that 
the light antiquark q will have zero momentum after scattering. It results in that the amplitude 
A will be divergent because the gluon propagator attached to q in Fig.l. To have meaningful 
predictions one can not neglect the space-time dependence here. The momentum of q is of order 
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of Aqcd, one can expand the amplitude in this momentum, this approach is used in [18-20]. We 
will use the approach here to regularize the divergence. For doing this we write the variables Xi 
as: 

xf^v- XiV^ + xf^i = LOiV^ + xf^i, for i = 1, 2, 3, (37) 

momenta can also be decomposed similarly as above. The dependence of the matrix element on 
XTiii — 1, 2, 3) can be safely neglected, and we have: 

' / ^ ^ ^ Tr.4ij (pi , p2 , A;i , ps) Tr (7°^^ (pi , P2 , A;2 , Pa) 7° ) ki 



with 



^1 = {rriQ + r)i)v, k2 = {niQ + r]2)v, h = rj^v. (39) 



To separate the divergences more clearly, we use translational invariance for the matrix element 
to shift the variable of fields with —UJ2V, and rearrange the variables as: 002 — > — <^2, a;3 — a;2 —>■ coi, 
UJ1—0J2 — > (^3, and r)i + Tjs — 7)2 ^ r)2, rji ^ Tjs, Tjs ^ rji. After this re-arrangement the momentum 
ki reads: 

h ^ {mq + r)3)v, A;2 = (mg + (773 + 771 - 772))^ , ks^rjiv. (40) 

The light quark q carries the momentum rjiv and r]2V in the amplitude A and respectively. 
Now we observe that for 773 the amplitude A and A^ is finite, while for t^i — > and 772 ^ 
they are divergent respectively. We expand the amplitude A and A^ respectively in 771^ and 7721' 
and only keep the leading terms in the expansion in 771, 772 and 773, i.e.; 

A, = ^(1 + ■■■), 4 = -(1 + ■■■), (41) 

771 772 

where T^j and T^j do not depends on 771, 772 and 773. Using Eq.(41) the integral over 773 and cos can 
be performed and we obtain: 



da 



(1) 



1 d^k f d^p. 



2s{27TfJ (27r)3' 

• J 2^ 977)772 """^"^'^ (Pi' P2, mgv, p3)Tr(7°T^ (pi , p2, rngv, P3)7°) w 

• J du^du2e+'^'^'-'^'-''^0\qk{u2v)hi{0)a^HQ^^^^ (42) 

The Dirac indices of the Fourier transformed matrix element can be decomposed. In general there 
are many terms. With the property in Eq.(33) the number of terms can be reduced greatly and 
we have only four terms: 
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- {Pv^tX m) - {i^iMu {PviT,i,\^ >vP(77i, m) (43) 



where 



^^^ji^-v -fv", = ^^^^ 



and the functions arc defined in the rest frame of Hq. By using matrix elements in the rest frame 
the factor in Eq.(43) appears because of the nonrelativistic normahzation of the state. The 
functions are: 

M'\m:V2) ^^e-'^-^^+'^-^^ ■ mu;,v)h{0)a^„^aH^h{0)q{u;,v)\0), 

M'\vi:V2) ^^e-'^-^^+'" ■ (0|g(a;2^)75/i(0)al,^a^^/.(0)75g(a;i^)|0), 

M'\m,m) = ^^6-"--+^^-- . (0|g>2^)7T/^(0)al,a^,/i(0)7T.?(a;i^)|0), 

M'\vuV2) = -1^/^^^"^'""^"'""' • (Ok>2^)7T75/i(0)al,^a^^/.(0)7T.75gK^)|0), 

(45) 

In the above definitions the operator qhq creates Hq in its rest frame. For the color-octet matrix 

element one has the same decomposition as in Eq.(43) with the functions W/^^ (771, 772, 773) with 
i = 1, ■ ■ ■ 4. The definition of these functions are obtained by replacing the color matrix 1 1 with 

rpa ^ rpa ^Yie above equation. 

To express our results as the differential cross-section like da/dt we need to consider the phase- 
space integral 

/(S-i§FP-)'*'&'^.+m-P.3). (46) 

Because Pi + P2 — ki — — ps we can insert the identity 

J dr]id{v ■ {pi +P2- mqv - r]iv - p^)) = 1, (47) 

where we neglect 773 in Eq.(40) and the momentum k of Hq is approximated by A; = mqv + rjiv. 
The (5-function in Eq.(47) combining those 5-functions in Eq.(46) gives usual 5-functions 5^(pi -|- 
P2 ~ k —793) for momentum conservation, then the phase-space integral can be calculated as usual. 
We obtain: 
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where we use k 



rnqv. Our results contain 8 parameters which represent nonperturbative effects. 

(^'«)(i = l,.--4): 



These parameters are defined as integrals with the functions W, 

(1,8) 1 [ dr]i dr]2 



w. 



rrVn J 



(49) 



These parameters are dimensionless. At first look, they may scale as Aq^j^/itlq. But, the dominant 
contribution of the integrals comes from the region where ?7i ~ 772 ~ Aqqd, refiecting the fact 
that the light antiquark q carries a small fraction of the momentum of Hq. This results in that 
the parameters wl^'^\i — 1, 2, 3, 4) scale as Aqcd/itiq. With these parameters our results for the 
process g + q ^ Hq + X are: 



where 



da 



(i) 
R 



dt 



(s,t} 



52 j i ' 



(50) 



(1) 



(1) 



-I 

81s t- 



127T2 - 16sT - QAf 



81s T 



64U 



81s 

1 + 



19 + 



r3 



C/2 



28T 



+ 



s 

368C/ 



Urp2 

64C/2 



T 



rp2 



2lP 
2^2 



48m^s(8s + 9T) . 



(51) 



(8) 
3 



bP 



2 rm^s(8s2 + 8sT + 9r2 



27s I- 
4s2 - sT + 4T2 



T3C/2 



2r 



79 



18T2 UU 

+ + 



[/2 



Bi'^ 



7717 



2T 



sU 

A; 

27s ^ 

233 + 



U 



[22 



9T 

266T2 316T 



18[/ 9T2 



[/2 



T 

lot/ 



f/2 

8t/2 



2^2 



2^2 



+ 



3m|s r 



T 



8U' 
2^2 



9T 



8U 

Y 



]}. 



(52) 



with 



m: 



T. 



(53) 



Similarly for the process 7 + — > i^Q + X, we have 



dt 



3=1 



(54) 



where 
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Bi = B2 = -rA^^SrTi^ + + ^ ^ ^ + + — - — - ^ iTT^j, 

4.S I- T3/^/2 \ „^ TP I T Q T3J ef/ J' 



1 rm^jSr^ ^ ^ 4(2 + fi;)s2 2(3 + 7k)s 

^3 + rri2 T-- ~'~ 



B3 = B4 



4s ^ t/2 



s3 4(2 + fi:)s2 2(3 + 7k)s Sk^T 



T3 r2 T s 

1 + ^) (1 + ^)]}' (55) 



with C*-^'^^ given in Eq.(25). In above results polarizations of initial states are averaged. 

Our results contain 8 parameters in general which are defined as integrals with matrix elements 
of HQET. They are universal, i.e., they do not depend on a specific process. With our results one 
can check that the ratio is a constant: 

-r^, t_^_iti . ~ const., (56) 



i.e., the ratio does not depend on s, hence the contribution from quark recombination is not 
suppressed by inverse of certain power of s, and it can give a significant contribution. This 
contribution is only suppressed by the parameters given in Eq.(49). 

The 8 parameters can be effectively reduced to 4 because of the property of Bi''^^ — Bs^'^-* and 
can show this property without explicit calculation. With the expansion in 
Eq.(41) the last diagram in Fig.l. will not contribute. Taking any contribution from the first 4 
diagrams in Fig.l, e.g., the interference between amplitudes from Fig. la and Fig. lb, the coefficient 



of Wi^^ can be written as a trace: 



Bi') ~Tr{---77-p2/---}, (57) 



where • • • denote polynomials of products of 7 matrices, I denotes 4x4 unit matrix coming from 
Eq.(43). Since we do not observe the polarization of initial hadrons, the spin of the initial light 
quark g(p2) is averaged, it results in the factor 7 • P2- The same leads to the coefficient of W2'^: 

^Tr{---(-75)7-P275---}, (58) 

it is clear that both coefficient are same. If initial hadrons are polarized, the coefficients are in 

general not the same. To illustrate this, we consider the case where only the hadron containing 
the light quark q{p2) is transversally polarized with the polarization vector Si^ and another initial 
hadron is unpolarized. In this case spin-dependent parts of perturbative coefficients, denoted as 
^si (^ — 1' 2, 3, 4), can be calculated by replacing 7 • p2 with icTf^^'y^pi^SI^, the contribution to the 
spin-dependent part of the cross section is a convolution with the so-called transversity distribution 
function instead of usual parton distribution functions [26]. With the replacement one can easily 
find that spin-dependent parts of perturbative coefficients have the properties B^i^^ = —B[^'^^ and 
'^s3^'* = —B^^^'^\ in contrast to the spin-independent part. Hence, if the spin-dependent parts of 
perturbative coefficients are not zero, we will have 8 nonperturbative parameters which lead to 
different contributions to the cross section. With symmetries one can also identify when these 
spin-dependent parts are zero in this case. With the rotational symmetry, the spin-dependent 
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part of differential cross-sections consists of two parts, one is proportional to St • k, while another 
is proportional to St • (P x k), where P is the three momentum of the polarized hadron. If the 
parity is conserved, the term proportional to St ■ k is zero. With the time-reversal symmetry one 
can show that the term proportional to St ■ (P x k) is zero if no absorptive parts in amplitudes. 
Further, one can show with arguments from chirality that the two terms are proportional to mq. 
Therefore, in the case considered here, we will have effectively 4 parameters only provided that 
there is no parity-violating interaction and absorptive part in partonic processes. Since parity- 
violating interactions can be involved and an absorptive part can also appear beyond tree-level, 
we have in general 8 parameters in differential cross-section. Similarly, one can also show that 
the 8 parameters lead to different contributions in the case if all initial hadrons are polarized. 
It should be noted that in the above discussion the polarization of the final hadron, if it is not 
spin-less, is summed. 

Because all these parameters are unknown, a numerical prediction can not be made. As an 
estimation, one may use vacuum saturation for wf'\i = 1, • • • ,4), i.e., the color singlet parts, to 
relate them to wave functions of Hq. However, the approximation of the vacuum saturation is not 
well established because g is a light quark, and also the approximation definitely does not apply 
for the color octet part, i.e., for wl^\i = 1, ■ • • ,4). With the definitions of these parameters one 
may use nonperturbative methods like QCD sum rules to study them, or they can be extracted 
from experiment. Once their numerical values are known, numerical predictions can be made with 
our results. In this work we do not make an attempt to fit experimental results with our results 
for determining these parameters, because in [18-20] it is shown that one can already describe the 
leading particle effect observed in experiment by keeping contributions with ^2 3 and neglecting 
other contributions. The relation between the parameters in [19] and ours can be identified: 



state, respectively. Our results of these contributions to the cross-section are exactly the same as 
those given in [19,20]. Prom our general analysis one can see that the pair Qq can form a scalar, 
a pseudoscalar, a vector, and a pseudovector state. All of them can lead to contributions to the 
leading particle effect, which can be predicted with the above results by using Eq.(28) or Eq.(29). 
If the light quark q is replaced with a heavy quark, one can use NRQCD factorization [21] to make 
prediction, then at leading order of the factorization, the Qq can only form a pseudoscalar or vector 
state, i.e., ^5*0- and ^y^i state. In this approach motivated by NRQCD factorization only terms 
with W2 3 will remain, while other terms are at higher order in the approach. This approach 
is used in [18-20]. With this approach for polarized beams one will still have 4 parameters in 
contrast to our results containing 8 parameters. 

4. Summary 

In this work we have made a general analysis for the effect of quark recombination in inclu- 




^Sn and ^S, 



(59) 



1 
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sive production of a heavy meson. Through these effects the leading particle effect observed in 
experiment can be explained. In the contributions of quark recombination to differential cross 
section we factorize perturbative- and nonperturbative parts in two different cases. One is for 
large transverse momentum of the produced hadron, where the transverse momentum is so large 
that any quark mass can be neglected. In this case we find the nonperturbative effect of quark 
recombination is parameterized by six functions defined with four quark operators at twist-4. 
Four of these functions contribute to leading particle effect, if the initial state is unpolarized. 
The observed effect of leading particle is of charmed hadron, one can take the charm quark as 
a heavy quark and perform the factorization with HQET. In this case the nonperturbative part 
consists of 8 parameters, which are defined as integrals of matrix elements with 4-quark operators 
in HQET. These parameters represent the transition of a quark pair Qq in different states into the 
produced hadron, the pair can be in color singlet or color octet state and forms different states as 
scalar, pseudoscalar, vector and pseudovector state. The results obtained here hold not only for 
large transverse momentum, but also for small transverse momentum. Our results are different 
than those obtained with a factorization motivated by NRQCD [18-20], in which the pair only 
forms pseudoscalar or vector state, although for unpolarized beams both results contain the same 
number of nonperturbative parameters. The perturbative part in both types of factorization are 
calculated at tree-level for hadro- and photoproduction. 

The nonperturbative functions or parameters are unknown yet. This fact prevents us from 
a numerical prediction for leading particle effect. But they can be studied by nonperturbative 
methods like QCD sum rules or models, or they can be extracted from experimental results. In 
this work we have not tried to extract the eight parameters from experimental results, because it 
is shown in [18-20] that one can already describe the leading particle effect observed in experiment 
only by keeping contributions of w^^^ , and contributions oi w{l in our result is proportional to 

the contributions of W2^^ , respectively. Hence, with effects of all possible states of Qq the leading 
particle effect can also be generated and can be described. It should be noted that the leading 
particle effect of charmed baryon can be also analyzed in a similar way presented in this work, 
results will be pubhshed elsewhere [27] . 
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